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ABSTRACT 

We present broadband X-ray spectral study of two ultraluminous X-ray sources (ULXs), M81 X-6 
and Holmberg IX X-1 based on Suzaku & XMM-Newton observations. We perform joint broadband 
spectral analysis of the brightest sources in the field, i.e. the two ULXs and the active galactic nucleus 
(AGN) in M81, and demonstrate that the X-ray spectra of the ULXs cut off at energies }^ 3 keV with 
negligible contribution at high energies in the Suzaku HXD/PIN band. The 90% upper limit on the 
10 — 30 keV band luminosity of an underlying broadband power-law component is 3.5 x 10^^ ergs s~^ 
for M81 X-6 and 1.2 x 10^® ergs s~^ for Holmberg IX X-1. These limits are more than an order of 
magnitude lower than the bolometric (0. 1 — 30 keV) luminosity of 6.8 x 10"^^ ergs s~^ for M81 X-6 and 
1.9 X 10''° ergs s^^ for Holmberg IX X-1. Our results confirm earlier indications of spectral cut-offs 
inferred from the XMM-Newton observations of bright ULXs and show that there isn't an additional 
high energy power-law component contributing significantly to the X-ray emission. The spectral form 
of the two ULXs are very different from those of Galactic black hole X-ray binaries (BHBs) or AGNs. 
This implies that the ULXs are neither simply scaled-up versions of stellar mass BHBs or scaled-down 
versions of AGNs. 

Subject headings: accretion, accretion disks — X-rays: binaries — X-rays: galaxies — X-rays: indi- 
vidual (M81 X-6, Holmberg IX X-1) 



1. INTRODUCTION 

Generally galaxies host numerous X-ray sources. Most 
off-nuclear, compact X-ray sources are believed to be 
X-ray binaries (XRBs) with a compact object - white 
dwarf, neutron star or a black hole. A few of the bright- 
est X-ray sources at the high end of the point source 
luminosity function can exceed the Eddington limit of 
even the most massive stellar mass black hole, sometimes 
by large factors. The nature of these sources, known as 
the ultra- luminous X-ray sources (ULXs), continues to 
remain an enigma since their dynamical mass measure- 
ments have not been possible. 

Possible explanations for the high luminosities of 
ULXs include {i) X-ray binaries with intermediate mass 
(Mhh ^ 10^ - 10^ Mrn) bla ck holes (IMBH) (e.g., 
IColbert fc Mushotzkvl lT999: Mil ler fc Colbe^l200l. (ii) 
X-ray binaries with anisotropic emission ( King et al.l 
|2001| ). {in) beamed XRBs with relativistic jets di- 
rectly point ing towards us i. e., scaled down versions 
of bl azars (jMirabel fc Rodriguez! [1999; Kording ct al. 
120021) , and (iv) XRBs with s uper-Eddington accretion 
rates (jBegelmanl I2001L 120021) . Chandra and XMM- 
Newton observations of ULXs have shown a variety of 
spectral shapes e.g., sim ple powerlaw similar to the 
low/hard state of BHBs ([Winter et all 120061 ) and soft 
excess (blackbody kT ~ 0.1 — 0.4 keV) plus power- 
law component. The soft X-ray excess emission has 
been interpreted as optically thick emission from a thin 
accretion disk with temperatures in the range of ~ 
100 — 300 eV, suggestive of an IMBH accreting at sub- 
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Eddington (^O.lL^rfrf) rates (iMil lcr et al. 2003. 2004b,^; 
iMiUer fc CoibertI [200l iFabian et al.. .2004 ). However, 
the cool disk plus powerlaw spectra of ULXs may 
not always cor respond to the high/soft state of BHBs 
([Roberts et all feOOS). High S/N XMM-Newton observa- 
tions of bright ULXs has revealed curva ture or cut-offs 
in th e spectra at hi gh energies ^ 5 keV (iStobbart et al.l 
200l iDewangan et al. 200i lAgrawal fc: Misral [20061: 
Gladstone fc Robertsi i2009i ) . However, it is not known 



if the X-ray spectrum actually has a cut-off or there is a 
deficit of emission near 7 keV due to perhaps an absorp- 
tion edge. Observation of strong X-ray emission above 
10 keV will put strong constraints on the nature of ULXs. 
The presence of a strong hard X-ray powerlaw will estab- 
lish ULXs to be either scaled-up BHBs or scaled-down 
versions of AGNs and the lack of emission near 7 keV 
may then be explained as the blurred iron K-edge due to 
Compton reflection arising from regions very close to a 
black hole. On the other hand, absence of such a strong 
hard X-ray component will make ULXs very different 
from BHBs or AGNs, probably accreting in a state that 
is not usually attained by them. 

Suzaku observations have shown that the broadband 
spect rum of M82 X-1 is slightly curved ([Mivawaki et al.l 
I2009D . However, M 82 X-1 is in a crowded field and it 
is not clear if the contribution of nearby multiple unre- 
solved sources both in the XIS and PIN FOV can al- 
ter the broadband spectral shape. Here we study the 
broadband spectral shape of two bright ULXs M81 X- 
6 and Holmberg IX X-1 (Ho IX X-1) located in the 
dwarf irregular galaxy Holmberg IX, a companion of 
M81. Ho IX X-1 is also known as M81 X-9. We con- 
strain the 0.6 — 30 keV spectral shape of the two ULXs 
for the first time based on joint spectral modeling of mul- 
tiple sources observed with Suzaku. 



2. OBSERVATIONS AND DATA REDUCTION 

We used the onaxis Suzaku observations of M81 (obser- 
vation ID: 906004010) and Holmberg IX X-1 (observa- 
tion ID; 707019010) performed during 15-16 September 
2011 and 13-17 April 2012 for exposure times of 45.6 ks 
and 182.5 ks, respectively. We also used XMM-Newton 
observation of Holmberg IX X-1 performed during 26-27 
September 2004 for an exposure time of 119.1 ks (obser- 
vation ID : 0200980101). 

We processed the unfiltcrcd Suzaku data with 
aepipeline available with HEASOFT version 6.13 and 
used the recent calibration database (CALDB) available 
as of 20 January 2013. We extracted the spectra of Ho 
IX X-1 from the filtered XIS event lists using a circu- 
lar region of radius 220", and the corresponding back- 
ground spectra from two circular regions each with ra- 
dius 110". The ULX M81 X-6 is 3.4' away from the 
AGN. Therefore, we extracted the XIS spectra for M81 
AGN and M81 X-6 using smaller circular regions of radii 
180" and 80", respectively. We extracted common back- 
ground spectra applicable to both AGN and ULX using 
two circular regions of radius 110" each. The circular 
extraction regions used to extract the spectral products 
are shown in Figure [TJ The XISO and XIS3 are the front 
illuminated (FI) CCDs and the two spectra for each ob- 
ject were combined using the tool addascaspec. We also 
extracted the HXD/PIN source and background spec- 
tra from both the observations using the task hxdpinxbpi 
from the PIN event lists and the tuned non X-ray back- 
ground files appropriate for the observations. 

Since the HXD/PIN is a non-imaging instrument with 
a field of view (FOV) of 34' x 34' FWHM, there may 
be contribution from sources within the FOV. To iden- 
tify the potential contaminating sources, we analyzed the 
XMM-Newton observation of Ho IX XOl. We processed 
the EPIC-pn and MOS data in a standard manner. We 
show the MOS image of the field containing Ho IX X-1 in 
Figure m In the MOS FOV of 30' diameter, apart from 
Ho IX X-1., there are two bright sources, M81 AGN and 
ULX M81 X-6 which can contribute to the PIN spec- 
trum. 

3. ANALYSIS & RESULTS 

We begin with the spectral analysis of the Suzaku ob- 
servation of M 81. While the AGN and the ULX are 
clearly resolved in the XIS data, the PIN data is addition- 
ally contaminated by the nearby bright ULX Ho IX X-1. 
Other sources seen in the XMM-Newton MOS image are 
too faint to contribute to the PIN spectrum. We fit- 
ted an absorbed powerlaw model to the PIN and the 
AGN XIS data. We used the the 16 - 30 keV band PIN 
data (as there is no detection above 30 keV) and the 
2 — 10 keV XIS band data whic h excludes the low energy 
thermal emission of the AGN (iPage et al.ll2003D. M 81 
AGN is known to show an iron line ( Dewangan et al.l 
|2004( ) and therefore we also used a Gaussian line. This 
simple model resulted in a good fit, x^ — 245 for 290 
degrees of freedom (dof ) . To find any contribution from 
the nearby bright ULXs to the HXD/PIN, we included 
an additional power-law component. We fixed the nor- 
malization of the second powerlaw component to zero for 
the XIS and varied it for the PIN data. This did not im- 
prove the fit and hence we calculated 90% upper-limit 



TABLE 1 
Best-fit parameters for the ULXs and AGN. 
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(") Model 1: AGN (wabsx (PL-fGL(l))) + ULX (wabsx cutoff pi); 
Model 2: AGN (wabsx (PL-fGL(l))) + ULX (wabs X (bbody 
+ cutoffpl + GL(2)); Model 3: AGN (wabsx (PL+GL(1))) 
+ ULX (wabsx (diskbb + nthcomp + GL(2)). C"' Galactic 

column in units of 10~-^0 cni~^. ^'^' Normalization in 

10"-^ photons cm~2 s"-' keV"-' at 1 keV. ''*' Line flux in 
10^* photons cm~2 s~^. ('^^ Normalization of cutoff PL in 
10^* photons cm~^ s~^ keV~^ at 1 keV. 

on the 10 — 30 keV flux by fixing the photon index r2 at 
different values ranging from 1 to 2.5. We find that the 
maximum possible contribution by the nearby ULXs is 
< 6.8 X 10"" ergs cm^^ g-i ^^ ^^le 10 - 30 keV band. 

To estimate the hard X-ray contribution from the 
ULXs, we fitted the XIS spectra of Ho IX X-1 and 
M81 X-6 with absorbed power-law models and calcu- 
lated the 10 — 30 keV band fluxes by extrapolating the 
best-fit models. We found /x(10 - 30 keV)= 8.3l^;2 x 
10^12 ergs cm^2 g-i foj. jjq jx X-1 and (3.6 ± 0.2) x 
10^12 foj. M81 X-6. Thus, the estimated fluxes for the 
ULXs based on the extension of the simple power-law 
models are much higher than the 90% upper-limit on 
the contribution by the nearby sources to the HXD/PIN 
data. This implies that the ULXs lack significant hard 
X-ray emission and their spectra cannot be described by 
a single power-law from soft-to-hard X-rays. This con- 
firms the earlier results of a high ene rgy cutoff in the 
XMM-Newton spectra of these ULXs (jDewangan et al.l 
[200l:[Stobbart et al.l 120061) . 

To better constrain the broadband X-ray spectral 
shapes of the ULXs, we performed joint spectral anal- 
ysis of the Suzaku XIS-I-PIN spectra of the ULXs and 
the AGN. For this purpose, we created ISIS fit functions 
which are a combination of models for the ULXs and the 
AGN but a given model component is applied to a spe- 
cific spectrum or a subset of spectra. First we fit the 
spectra obtained from the onaxis Suzaku observation of 
the M81 AGN. We used the combined XIS FI spectra 
of the M81 ULX and the AGN and the corresponding 
HXD/PIN spectrum of the field. We did not use XIS 
spectrum of Ho IX X-1 as it is not in the field of view. 
We also did not account for its contribution to the PN 
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Fig. 1. — Suzaku XIS images of Holmberg IX X-1 (left panel) and M81 
of the region containing Ho IX X-1, the ULX and AGN in M81. Also 
dashed hne for background that were used for spectral extraction. 

spectrum. This ULX is 12.4' away from the AGN in 
M81. Thus, given the triangular response of HXD/PIN 
with an FWHM of 34', the ULX contributes - 50% of 
its actual flux to the PIN data. As inferred above, Ho 
IX X-1 generally lacks the hard X-ray emission and its 
contribution to the PIN band is likely to be very small. 
We further deal with this issue below when we analyze 
the on- axis Suzaku observation of Ho IX X-1. 

For the joint spectral analysis of the M81 ULX and 
the AGN spectra, we used a joint model which is a com- 
bination of the model for the ULX (wabsx (cutoff pi 
+ gauss)) and the AGN (wabsx (powerlaw + gauss)). 
We refer to this as model 1. The ULX model only is ap- 
plied to XIS spectrum of M81 X-6 and the AGN model 
alone applied to the XIS spectrum of the AGN but both 
models are combined to flt to the PIN data. As before, 
we excluded XIS data below 2 keV for the AGN only 
and used a constant multiplicative factor of 1.16 for the 
PIN data. The model provided a statistically accept- 
able fit with x^ — 418.4/435 and the best-flt parameters 
are listed in in Table [T] while the unfolded spectral data 
and the best-fit model are shown in Figure [H We also 
examined the continuation of the powerlaw in the hard 
band by fixing Ecutoff at 100 keV. The fit became worse 
with x^ /dof = 679/436, ruling out the presence of such 
a powerlaw. Further, we added an additional powerlaw 
component with F = 2 to the best-fit ULX model 1 and 
estimated the 90% upper-limit on the 10 — 30 keV flux to 
be fx = 2.3 X 10~^^ ergs cm~^ s~^ and the correspond- 
ing luminosity Lx = 3.5 x 10^^ ergs s""'^ for an underlying 
power-law component. We also tested physical models 
that do not contribute signiflcantly in the hard X-ray 
band. We replaced the cutof fpl in model 1 with the 
standard disk model diskbb which resulted in a good 



AGN and ULX (middle panel) and the XMM-Newton MOS image 
shown are the circular regions with solid lines for the sources and 



1.62 



-1-0.04 
-0.05 



(3.7 ± 0.1) X 10-12 ergs cm'^ s 



keV, 

1 



fit ix^/dof = 432.4/436) with kTi 

/x(0.1 - 10 keV) 
We also tested the p-free disk model diskpbb which al 
lows the disk temperature to scale with radius as r~P. 
For standard disk model, p = 0.75 an d lower values 
of p may indicate a slim accretion disk (jWatarai et al.l 



l2000f ). The p-free disk model resulted in a slightly 
better fit {x^/dof = 421.5/435) with kT„, = 1.8 ± 
0.1, p = 0.69 ± 0.02 and /x(0.1 - 30 keV) = (4.3 ± 
0.1) X 10-12 ergs cm^^ g-i and ix((0.1 - 30 keV) = 



6.8 X 10^^ ergs s^^. The parameters of the AGN model 
components remained similar for different ULX models 
cutof fpl, diskbb and the p-free disk. Thus, we confirm 
that the M81 X-6 spectrum cuts-off at high energies and 
does not contribute significantly to the HXD/PIN band. 
Moreover, any possible hard X-ray contribution by the 
ULX Ho IX X-1 (which we show below is unlikely) would 
further strengthen the result that M81 X-6 is dim in hard 
X-rays. 

We next perform the joint spectral analysis of the on- 
axis Suzaku observation of Ho IX X-1. The ULX is 
well isolated spatially and its X-ray emission in the XIS 
band is not contaminated by any other strong source 
(see Fig. [Ij. However, the PIN data is obviously con- 
taminated by the AGN in M81. We used a joint spec- 
tral model similar to the model 1 that combines the 
ULX model (wabsx (bbody + cutof fpl + gauss) and 
the AGN model (wabsx (powerlaw + gauss)). We re- 
fer this as model 2. As before we used a relative normal- 
ization of 1.16 between the XIS spectrum of Ho IX X-1 
and its contribution to the PIN spectrum. Since there is 
no simultaneous XIS spectrum of the AGN, we used the 
AGN XIS spectrum from the onaxis observation of M81. 
We used a variable relative normalization between the 
XIS spectrum of the AGN and its contribution to the PIN 
band to allow for flux variability of the AGN between the 
two Suzaku observation. Thus, the use of the XIS AGN 
spectrum in modeling the ULX spectrum is equivalent 
to assuming that the AGN did not change its spectral 
shape. This is justifled as the 2 — 10 keV spectral shape 
of M81 AGN is known t o show very little variability if at 
all (jlshisaki et al.|[T996() . As before we excluded the XIS 
spectrum of AGN below 2 keV. The joint fit resulted 
in a statistically acceptable fit {x^/dof — 509.3/517). 
The best-fit parameters are listed in Table [1] while The 
unfolded spectrum and the best-flt models for the ULX 
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Fig. 2. — Unfolded XIS and PIN spectra of ULXs and AGN. Lejt panel: The XIS and PIN spectral data and the best-fit models for the 
AGN and the ULX in M81 derived from the onaxis Suzaku observation of M81 on 15 September 2011. Right panel: The XIS spectrum 
of M81 AGN obtained from the Suzaku observation of 15 September 2011, the XIS spectrum of Ho IX X-1 and the corresponding PIN 
spectrum derived from the onaxis Suzaku observation of Ho IX X-1 on 13 April 2012. Also shown are the individual spectral component — 
blackbody, Gaussian fine and the cutoffpl for Ho IX X-1. 

we showed that the broadband spectra of M81 X-6 and 
Ho IX X~l cutoff at 2.8 keV and 8 keV, respectively. 
This conforms the high energy spectral curvature in these 
ULXs ea rlier found with the limited bandpass of XM M- 
Newton ()Dewangan et ani2006l: IStobbart et al.ll2006D . In 
addition, we do not find additional high energy spec- 
tral components similar to that observed from BHBs and 
AGNs. This demonstrates clearly that the two ULXs 
are in unusual spectral states that are distinctly differ- 
ent than that observed from BHBs and AGNs. 

The spectrum of M81 X-6 is well described by the 
standard disk model with high temperature fcT,:„ = 
1.621q'q5 keV similar to that observed from BHBs but 



and AGN are shown in Figure [21 The best-fit parame- 
ters for the AGN were similar to that obtained earlier 
(see Table [1]). As expected from the earlier analysis, the 
spectrum of Ho IX X-1 cuts-off and makes very little 
contribution to the HXD/PIN band. A weak iron line 
is detected for Ax^ = —11.5 for two additional parame- 
ters i.e., at a significance level of 99.7% according to the 
F-test. The line width was not well constrained, there- 
fore a was fixed at its best-fit value. Such a weak iron 
line wa s not detected in the long XMM-Newton obser- 
vation (jDewangan et al.ll2006f) and is most likely associ- 
ated with the host galaxy. The relative normalization 
for the AGN contribution to the PIN data is CSSlg-jg. 
Fixing the value to 1.16 does not change the cutoff en- 
ergy for the ULX and hence uncertainties in the AGN 
flux does not effect the overall result. We further tested 
the reality of the spectral cut-off of Ho IX X-1 by fixing 
the cutoff energy at 100 keV. This resulted in a poor 
fit (Ax^ — -1-99.4 for the loss of only one parameter). 
Thus, the presence of spectral cut-off is highly signif- 
icant. As before, we estimated 90% upper limit on the 
10 — 30 keV flux of an underlying powerlaw component to 
be fx = 7.5 X 10"-^'^ ergs cm~^ s~^ and the correspond- 
ing luminosity to be Lx = 1-2 x 10'^^ ergs s~^. We also 
tested the thermal Comptonization model by replacing 
the cutoffpl in model 2 with the nthcomp and bbody 
with the diskbb component (model 3), the joint fit re- 
sulted in x^/dof = 519/517) with kT^ = 2.6 ± 0.2 keV. 
In this fit, we treated the diskbb component with kTin ~ 
0.35 keV to be the seed photon source for the Comp- 
tonization. We calculated the flux of the diskbb and 
nthcomp components to be 5.4 x 10~^^ ergs cm~^ s~^ 
and 1.2 x 10~^^ ergs cm~^ s~^, respectively. 

4. DISCUSSION 

Using the broadband Suzaku observations of the fields 
containing two bright ULXs, we measured X-ray lumi 
nosities of 6.8 x lO^^ ergs s^^ for M81 X-6 and 1.9 x 10" 
for Ho IX X-1. With the joint spectral analysis of XIS 
and PIN spectra of the brightest sources in the field. 



the inferred X-ray luminosity exceeds the Eddington lu- 
minosity of AAMq black hole. Generally BHBs show 
standard disk spectra in their high/soft state with the rel- 
ative accretion rate in the range of ~ 0.1 — 0.5 ()Esin et al.l 
I1997D . Thus, the observed X-ray luminosity of M81 X-6 
will correspond to ^ IOOMq. The temperature of the 
standard disk around such a large black hole is unlikely 
to exceed 1 keV. Thus, M81 X-6 is not in the high/soft 
state generally observed in BHBs. The p-free disk fit 



showed that the temperature of disk [T oc 



-0.68 



) at dif- 



ferent radii is lower than that expected from a standard 



disk [T oc 



,-0.75 



). This suggests that the ULX M81 X-6 
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likely hosts a slim accretion disk. 

We confirmed the soft excess from Ho IX X-1 and 
showed that the spectral downturn near 8 keV earlier 
observed with XMM-Newton continues in the hard X- 
ray band above 10 keV. The soft excess and the spectral 
turnover is well described by a standa rd cool disk and 
cool Comptonization . as fo und earlier (jDewangan et aLl 
120061: IStobbart et all [20061 ). Based on the spectral anal- 
ysis of high quality X-ray spectra of ULXs observed with 
XMM-Newton, Gladstone & Roberts (2009) showed that 
most ULX spectra are best described by an outer cool 
disk and optically thick, cool Comptonization. They 
termed such a spectral form as the ultraluminous state 
which is thought to represent spectral transition between 
Eddington and super-Eddington accretion flows. It is in- 
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teresting to note that the luminosity of the soft excess 
component, described as diskbb, is 8.4 x 10"^* ergs s~^ 
corresponding to relative accretion of m/m'edd ^ 0.2 for 
a SOMq black hole. Thus, if there is mild beaming of the 
Comptonized component, a stellar mass black hole can 
explain the observed luminosity of Ho IX X-1. 

Finally, we note that while many BHBs and 
AGNs show strong relativistically broadened iron line, 
we are yet to observe such l ines from any ULX. 
iCaballero-Garcia fc FabianI (|201Cl[ ) explained the soft ex- 
cess and drop near 7 keV in the spectra of several ULX 



including M81 X-6 and Ho IX X-1 in terms of relativis- 
tically blurred Compton reflection from accretion disks. 
Such a model requires strong continuum above 7.1 keV. 
However, our finding of the spectral downturn and its 
continuation in the hard band above 10 keV in Ho X-1 
and M81 X-6 imply the absence of a strong illuminating 
continuum. Thus, our findings rule out the possibility of 
explaining the spectral downturn near 7 keV as a blurred 
iron K-edge. The absence of a strong illuminating contin- 
uum is most likely the reason behind the non-detection 
of relativistically broadened iron lines from ULXs. 
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